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ABSTRACT 

We present an analysis of the data produced by the MaNGA prototype run (P-MaNGA), aiming 
to test how the radial gradients in recent star formation histories, as indicated by the 4000A-break 
(D„(4000)), H<5 absorption (EW(Hd J 4 )) and Ha emission (EW(Ha)) indices, can be useful for under¬ 
standing disk growth and star formation cessation in local galaxies. We classify 12 galaxies observed 
on two P-MaNGA plates as either centrally quiescent (CQ) or centrally star-forming (CSF), according 
to whether D„(4000) measured in the central spaxel of each datacube exceeds 1.6. For each spaxel we 
generate both 2D maps and radial profiles of D„(4000), EW(Hd y i) and EW(Ha). We find that CSF 
galaxies generally show very weak or no radial variation in these diagnostics. In contrast, CQ galaxies 
present significant radial gradients, in the sense that D„(4000) decreases, while both EW(H^) and 
EW(Ha) increase from the galactic center outward. The outer regions of the galaxies show greater 
scatter on diagrams relating the three parameters than their central parts. In particular, the clear 
separation between centrally-measured quiescent and star-forming galaxies in these diagnostic planes 
is largely filled in by the outer parts of galaxies whose global colors place them in the green valley, 
supporting the idea that the green valley represents a transition between blue-cloud and red-sequence 
phases, at least in our small sample. These results are consistent with a picture in which the cessation 
of star formation propagates from the center of a galaxy outwards as it moves to the red sequence. 

Subject headings: galaxies:general - galaxies:stellar content - galaxies:formation - galaxies:evolution 
- surveys:galaxies - methods:observational 
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1. INTRODUCTION 

Our understanding of galaxies in the local universe has 
improved dramatically over the past decade, in signifi¬ 
cant part due to large optical spectroscopic surveys such 
as the Two-degree Field Galaxy Redshift Survey (2dF- 
GRS; Colless et al. 2001) and the Sloan Digital Sky Sur¬ 
vey (SDSS; York et al. 2000). One of the major find¬ 
ings of the SDSS is the discovery of ‘galaxy bimodal¬ 
ity’, through which local galaxies are divided into two 
distinct populations, termed ‘red sequence’ and ‘blue 
cloud’, according to their rest-frame colors or specific star 
formation rate (e.g. Strateva et al. 2001; Blanton et al. 
2003; Kauffmann et al. 2003a; Baldry et al. 2004). Red- 
sequence galaxies typically present bulge-dominated 
morphology, red optical/UV color, and little gas and star 
formation, while blue-cloud galaxies are usually gas-rich 
with blue colors, ongoing star formation and disky mor¬ 
phology. The smaller number of galaxies found lying in 
the ‘green valley’ between these two main populations in 
color space are widely believed to be caught in a tran¬ 
sition phase from the blue cloud to the red sequence 
(e.g. Bell et al. 2004; Faber et al. 2007; Martin et al. 
2007; Schawinski et al. 2007; Schiminovich et al. 2007; 
Wyder et al. 2007; Mendez et al. 2011; Gongalves et al. 

2012) , although this picture has recently been dis¬ 
puted by some authors (e.g. Schawinski et al. 2014). 
Narrow-field deep surveys reveal that this color bi¬ 
modality persists out to redshift of at least 2.5 
(e.g. Bell et al. 2004; Bundy et al. 2006; Cirasuolo et al. 
2007; Cooper et al. 2008; Faber et al. 2007; Martin et al. 
2007; Cowie & Barger 2008; Brammer et al. 2009; 
Williams et al. 2009; Muzzin et al. 2012; Huang et al. 
2013a), though the prevalence of the red population has 
significantly increased since a redshift of unity (Bell et al. 
2004; Bundy et al. 2006; Faber et al. 2007). The cessa¬ 
tion of star formation is thus an important process which 
has been contributing to galaxy evolution over the past 
~ 8 Gyr. 

A complete picture of the way in which the star 
formation in galaxies gets shut down remains elusive. 
However, recent studies of the scaling relations of galaxy 
properties and their dependence on local environment 
have clearly established that, in addition to stellar 
mass, both internal structural properties and external 
environment are key indicators, or may even be drivers, 
of the star formation cessation processes in galaxies 
(e.g. Kauffmann et al. 2006; Bell 2008; Franx et al. 
2008; Peng et al. 2010; Thomas et al. 2010; Bell et al. 
2012; Cheung et al. 2012; Fabello et al. 2012; Li et al. 
2012b; Fang et al. 2013; Mendel et al. 2013; Zhang et al. 

2013) . For instance, when studying the central galaxies 
in groups or clusters, the presence of a prominent 
bulge-like structure is found to be a necessary (but not 
sufficient) condition for stopping star formation (Bell 
2008; Bell et al. 2012; Cheung et al. 2012; Fang et al. 
2013). In addition, studies of the relationship be¬ 
tween galaxy morphology and color have revealed a 
significant population of red-sequence galaxies with 
disc-dominated spirals at both low-z (e.g. Wolf et al. 
2005, 2009; Bamford et al. 2009; Masters et al. 2010) 
and 2 ~ 1 — 2 (e.g. Bundy et al. 2010), which are 
preferentially found in galaxies with large bulges 
(Bundy et al. 2010; Masters et al. 2010). These findings 


support the ‘morphological quenching’ mechanism 
proposed by Martig et al. (2009), although studies of 
cold gas in massive galaxies indicate that a reduction 
in gas content is also required (e.g. Fabello et al. 2011). 
For central galaxies in massive dark matter halos 
above a critical mass of ~ 10 12 Mq, “radio-mode” 

AGN feedback and shock-heating may effectively re¬ 
duce the gas cooling efficiency, thus also playing an 
important role in preventing further star formation 
(Silk 1977; Rees & Ostriker 1977; Blumenthal et al. 
1984; Birnboim & Dekel 2003; Keres et al. 2005; 
Dekel & Birnboim 2006; Cattaneo et al. 2006). Pow¬ 
erful AGN feedback can also happen in the so-called 
“quasar mode”, which is predicted to be triggered by 
major mergers of gas-rich galaxies with comparable mass 
(e.g. Di Matteo et al. 2005; Hopkins et al. 2006), and is 
commonly adopted as one of the quenching processes in 
semi-analytic models. 

For satellite galaxies, star formation quenching seems 
to be driven primarily by external processes occur¬ 
ring within their host group/cluster, such as gas strip¬ 
ping by ram-pressure (Gunn & Gott 1972; Abadi et al. 
1999) and tidal interactions (Toomre & Toomre 1972; 
Moore et al. 1996). The ‘smoking gun’ of observations 
showing Hi gas being removed by ram-pressure strip¬ 
ping have been obtained for spirals in nearby clus¬ 
ters of galaxies (e.g. Vogt et al. 2004; Chung et al. 2009; 
Merluzzi et al. 2013). Studies of color profiles and sur¬ 
face brightness profiles for satellite galaxies in SDSS 
group systems prefer more gentle processes such as 
“starvation” or “strangulation” (e.g. Weinmann et al. 
2009). On the other hand, a recent study of the spe¬ 
cific star formation rate-stellar mass relation for galaxy 
groups/clusters at higher redshift (0.2 < z < 0.8, 
Lin et al. 2014) suggests that galaxy mergers play a pri¬ 
mary role in quenching satellites in galaxy groups, while 
strangulation is a process more important in cluster-scale 
environment. There have also been statistical studies 
of the cold gas content and star formation in galaxies 
residing in different environments, using either direct 
atomic gas mass measurements from large blind surveys 
of Hi emission such as the Hi Parkes All Sky Survey 
(HIPASS; Zwaan et al. 2005) and the Arecibo Legacy 
Fast ALFA (ALFALFA; Giovanelli et al. 2005) survey, 
or indirect estimates of gas content from the optical 
photometry in surveys like SDSS (e.g. Li et al. 2012b; 
Zhang et al. 2013). These studies have revealed that gas- 
related quenching depends not only on the stellar mass 
of the satellite galaxies, but also on their surface mass 
density. 

A full understanding of star formation cessation in 
galaxies thus requires spatially-resolved measurements 
of stellar and gaseous components to be obtained for 
a large sample of galaxies covering wide ranges in 
mass and color, probing a range of environmental 
conditions. Previous studies for this purpose have 
been mainly based on multi-wavelength broadband pho¬ 
tometry, for both nearby galaxies (e.g. de Jong 1996; 
Bell & de Jong 2000; Taylor 2005; Munoz-Mateos et al. 
2007; Zibetti et al. 2009; Roche et al. 2010; Suh et al. 
2010; Tortora et al. 2010; Gonzalez-Perez et al. 2011; 
Tortora et al. 2011; Kauffmann 2015) and those at high 
redshifts (e.g. Abraham et al. 1999; Azzollini et al. 2009; 
Szomoru et al. 2012; Wuyts et al. 2012; Szomoru et al. 
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2013; Wuyts et al. 2013), while some authors also made 
use of long-slit spectroscopy (e.g. Moran et al. 2012; 
Huang et al. 2013b). Significant improvements have 
been made in the past decade, thanks to many in¬ 
tegral field unit (IFU) surveys which have obtained 
spatially resolved spectroscopy for samples of galax¬ 
ies. However, to-date, the samples have been rela¬ 
tively small (< 100 galaxies): SAURON (Bacon et al. 

2001), DiskMass (Bershady et al. 2010), ATLAS'® 
(Cappellari et al. 2011). 

The ongoing CALIFA survey (Calar Alto Large In¬ 
tegral Field Area Sanchez et al. 2012a) is observing a 
sample of 600 galaxies, making a big step forward 
in 2D spectroscopic studies of nearby galaxies. For 
example, this survey has, for the first time, quanti¬ 
fied the spatially-resolved history of stellar mass as¬ 
sembly for galaxies beyond the Local Group, demon¬ 
strating how massive galaxies grow their stellar mass 
inside-out (Perez et al. 2013). Such an inside-out pic¬ 
ture of galaxy formation is supported by further mea¬ 
surements of gradients in oxygen abundance of Hu re¬ 
gions (Sanchez et al. 2014), in stellar metallicity of face- 
on spirals (Sanchez-Blazquez et al. 2014), and in stellar 
age and local surface mass density of galaxies of differ¬ 
ent morphologies and masses (Gonzalez Delgado et al. 
2014), all from the CALIFA survey. However, the sam¬ 
ple is still relatively modest, limiting the range of en¬ 
vironments and galaxy properties that can be explored. 
The more recently launched SAMI survey (Croom et al. 
2012) takes a further step with a goal of obtaining IFU 
observations for 3400 galaxies by 2016. In the meantime, 
the KMOS 3Z:i Survey is pushing forward the IFU obser¬ 
vations for high redshift galaxies, observing a sample of 
600 galaxies at 0.7 < z < 2.7 using KMOS at the Very 
Large Telescope (Wisnioski et al. 2015). 

As a next step at low-z, the study of these phe¬ 
nomena is a key goal of the upcoming MaNGA (Map¬ 
ping Nearby Galaxies at Apache Point Observatory, 
Bundy et al. 2015) survey. As one of the major pro¬ 
grams of the fourth-generation of SDSS, MaNGA will 
obtain two-dimensional, integral-field spectroscopy (IFS) 
for 10,000 galaxies in the nearby redshift range 0.01 < 
2 < 0.15, optimally selected for uniform coverage and 
resolution using the single (central) spectra from SDSS. 
Each galaxy will be spectroscopically observed with an 
IFU to obtain multiple high signal-to-noise ratio (S/N) 
spectra across the full optical wavelength range from 
3600 to 10300A. MaNGA will thus provide 2D maps of 
stellar populations and recent star formation histories for 
a large number of quiescent and star-forming galaxies, 
thus enabling extensive exploration of the various pro¬ 
cesses by which star formation ceases in the local uni¬ 
verse. 

In this paper, to assess the potential use of the MaNGA 
survey in studying the cessation of star formation, and 
to derive some initial results in this area, we present an 
analysis of the MaNGA prototype (P-MaNGA) data for 
a set of 12 galaxies, which were observed in the pilot run 
in 2012 December and 2013 January. These observations 
were made through time granted by the Sloan Digital Sky 
Survey-Ill (SDSS-III, Eisenstein et al. 2011). The proto¬ 
type data are described in more detail in Bundy et al. 
(2015), which also provides a general overview of the 


MaNGA project. 

For each of the P-MaNGA galaxies and for each spec¬ 
trum in the associated data cube, we have performed 
full spectral fitting using the methodology developed in 
Li et al. (2005), in order to separate out emission line 
components from the underlying absorption-line stellar 
spectrum. From these decomposed spectra, we have ob¬ 
tained full 2D maps of three key parameters- D„(4000) 
(the depth of the 4000A break), EW(HJa) (the equiva¬ 
lent width of the H<5 absorption line) and Ha (the equiv¬ 
alent width of the Ha emission line)—which are known 
to be sensitive indicators of stellar populations of differ¬ 
ent ages (Kauffmann et al. 2003b; Kauffmann 2014). We 
then investigate any radial variations in the recent star 
formation histories of the galaxies by plotting the val¬ 
ues of these parameters for different subregions of each 
galaxy, to see how they change with position. As we will 
show, the P-MaNGA sample, though small in size, re¬ 
veals intriguing systematic trends in the three diagnostic 
parameters across these galaxies, demonstrating the po¬ 
tential for using MaNGA (and similar IFU surveys) to 
understand disk growth and star formation cessation in 
local galaxies. 

Throughout this paper we assume a A cold dark mat¬ 
ter cosmology model with f l m = 0.27, Da = 0.73 and 
h = 0.7. A Chabrier (2003) stellar initial mass function 
(IMF) is adopted. 

2. DATA 

2.1. P-MaNGA Observations 

In the present work, we use data obtained using the 
MaNGA engineering prototype instrument (hereafter P- 
MaNGA) in 2013 January. P-MaNGA was designed to 
explore a variety of instrument design options, observ¬ 
ing strategies, and data processing algorithms. These P- 
MaNGA data offer valuable insights into MaNGA’s po¬ 
tential, but differ substantially from the MaNGA survey 
data in several ways. First, P-MaNGA used just one of 
the two BOSS spectrographs, with only 560 total fibers 
and 8 IFUs with three sizes: 19 fibers (Nig), 61 fibers 
(N 6 i), and 127 fibers (N 12 r). The P-MaNGA IFU com¬ 
plement was dramatically different from the MaNGA sur¬ 
vey instrument, with 5xNig (instead of just two), lxNgi 
(instead of 4), 2xNi 2 7 (instead of five), and no 37-fiber 
or 91-fiber IFUs (see Drory et al. 2015, for detailed de¬ 
scription of MaNGA instrumentation). Finally, we also 
note that the masses and sizes of the P-MaNGA target 
galaxies were not selected to be representative of the full 
MaNGA sample. 

P-MaNGA observations were obtained of three galaxy 
fields (Field 9, 11, 4) using SDSS plates 6650, 6651, 
and 6652, each containing 6 galaxies. Some of the P- 
MaNGA targets were drawn from early versions of the 
full MaNGA sample (Wake et al., in prep), but in many 
cases P-MaNGA targets were chosen for specific reasons. 
In each of the three plates, one N127 IFU was allocated for 
comparison purposes to a galaxy observed by the CAL¬ 
IFA survey (Sanchez et al. 2012b), even if the galaxy 
would not otherwise satisfy the MaNGA selection cuts. 
Additionally, the non-optimal IFU complement of the P- 
MaNGA instrument required some targets to be selected 
manually. For each plate, observations were obtained in 
sets of three 20-minute exposures, which were dithered by 
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TABLE 1 

Prototype Run Galaxies Observed 


BundleJField P-MaNGA 
Name 

mangalD 

IAU Name 

z 

log 

M* 

Mi 

mag 

(fl - r) 
mag 

Re 

typef 

7?ifu 

(Re) 




Field 9: PlatelD = 6650, 

3.0 hr, seeing V! 

7 





ma003_9 

p9-127A 

8-188794 

J093457.30+214220.9 

0.013 

10.7 

-21.3 

0.70 

23.7 

-k 

0.7 

ma008_9 

p9-127B 

8-131835 

J093506.31+213739.5 

0.013 

9.1 

-18.5 

0.51 

6.8 

i 

2.4 

ma002_9 

p9-61A 

8-188807 

J093712.30+214005.0 

0.019 

10.1 

-20.5 

1.30 

9.3 

•k 

1.2 

ma005_9 

p9-19B 

8-131893 

J093109.60+224447.4 

0.051 

10.6 

-22.2 

0.86 

4.0 

i 

1.6 

ma007_9 

p9-19D 

8-131577 

J093109.07+205500.5 

0.034 

10.3 

-21.2 

0.79 

3.2 

i 

2.0 

ma001_9 

p9-19E 

8-131821 

J094030.22+211513.7 

0.024 

9.7 

-20.0 

0.73 

2.9 

i 

2.2 



Field 4: PlatelD = 6652, 

2.0 hr, seeing 1'! 

3 





ma003_4 

p4-127A 

8-109661 

J105555.26+365141.4 

0.022 

10.7 

-22.1 

0.84 

10.3 

★ 

1.6 

ma008_4 

p4-127B 

8-109682 

J105259.05+373648.2 

0.042 

11.0 

-22.7 

0.77 

13.9 

★ 

1.2 

ma002_4 

p4-61A 

8-113576 

J105746.61+361657.8 

0.030 

9.7 

-20.0 

0.71 

3.4 

~k 

3.4 

ma004_4 

p4-19A 

8-113557 

J110012.10+362313.8 

0.027 

9.4 

-19.5 

0.49 

2.5 

i 

2.5 

ma005_4 

p4-19B 

8-113506 

J104958.69+362454.0 

0.023 

9.5 

-19.5 

0.53 

4.6 

i 

1.4 

ma006_4 

p4-19C 

8-109657 

J105605.68+365736.1 

0.022 

9.5 

-19.4 

0.72 

4.8 

i 

1.3 


Note. — The P-MaNGA Name is composed of a “p” for “prototype” followed by the field ID, a hyphen, and then 
the shorthand ID for the bundle used. This shorthand includes a number corresponding to Njpy for that bundle. The 
(g — r) magnitude is extinction-corrected. 

I Target type 1 indicates the galaxy would be selected in the MaNGA Survey’s Primary Sample. These galaxies often 
have spatial coverage to larger than 1.5f? e due to the different bundle size distribution used in P-MaNGA. A indicates 
a galaxy that was chosen by hand for the prototype run and does not pass the nominal selection criteria. 


approximately a fiber radius along the vertices of an equi¬ 
lateral triangle to provide uniform coverage across each 
IFU. These three plates were observed to varying depths 
and in varying conditions, as required by the P-MaNGA 
engineering tasks. Although plate 6650 (Field 9) was ob¬ 
served to a depth comparable to what will be regularly 
achieved during MaNGA operations, plates 6651 (Field 
11) and 6652 (Field 4) are both significantly shallower 
than MaNGA survey data, and plate 6651 was inten¬ 
tionally observed at high airmass (~ 1.5) resulting in 
particularly poor image quality. In this paper we use 
data obtained with plates 6650 and 6652 (see Table 1), 
while excluding plate 6651. 

The raw data were reduced using a prototype of the 
MaNGA Data Reduction Pipeline (DRP), which is de¬ 
scribed in detail by Law et al. (2015). In brief, individ¬ 
ual fiber flux and inverse variance spectra were extracted 
using a row-by-row algorithm, wavelength calibration 
was obtained from a series of Neon-Mercury-Cadmium 
arc lines, and the spectra were flatfielded using internal 
quartz calibration lamps. Sky-subtraction of the IFU 
fiber spectra was performed by constructing a cubic ba¬ 
sis spline model of the sky background flux as seen by the 
41 individual fibers placed on blank regions of sky, and 
subtracting the resulting composite spectrum shifted to 
the wavelength solution of each IFU fiber. 

Flux calibration of the P-MaNGA data is performed 
by fitting Kurucz model stellar spectra to the spectra of 
calibration standard stars covered with single fibers at 
each of the three dither positions. The flux calibration 
vectors derived from these single-fiber spectra were found 
to vary by ~ 10% from exposure to exposure, depending 
on the amount of light lost from the fiber due to at¬ 
mospheric seeing and astrometric misalignments. While 
this uncertainty is acceptable for the present science pur¬ 
poses, the flux calibration uncertainty of the single fibers 
ultimately drove the design decision of the MaNGA sur¬ 


vey to instead use 7-fiber IFU “mini-bundles” for each 
calibration standard star, which results in a few percent 
photometric uncertainty (see Yan et al. 2015). 

Flux-calibrated spectra from the blue and red cam¬ 
eras were combined across the dichroic break using an 
inverse-variance weighted basis spline function. Astro- 
metric solutions were derived for each individual fiber 
spectrum that incorporate information about the IFU 
bundle metrology (i.e., fiber location within an IFU), 
dithering, and atmospheric chromatic differential refrac¬ 
tion, among other effects. Fiber spectra from all expo¬ 
sures for a given galaxy were then combined into a single 
datacube (and corresponding inverse variance array) us¬ 
ing these astrometric solutions and a nearest-neighbor 
sampling algorithm similar to that used by the CALIFA 
survey. For the P-MaNGA datacubes, a spaxel size of O'.'S 
was chosen. The typical effective spatial resolution in the 
reconstructed datacubes can be described by a Gaussian 
with a full width at half maximum FWHM ss 2"5. When 
binning the datacubes, we scale the resulting error vec¬ 
tors to account, at least approximately, for wavelength 
and spatial covariance in the P-MaNGA error cubes. 

2.2. SDSS-I perspective of the P-MaNGA galaxies 

All our galaxies have photometry in the SDSS-I five 
bands, as well as archival spectroscopy from the SDSS-I 
3'/0-diameter fiber. Figures 1 and 2 show the properties 
of this sample based on the SDSS data. With these two 
figures we aim to show how our galaxies compare to the 
larger SDSS-I sample from which they are drawn. 

Figure 1 presents the galaxies on the plane of stel¬ 
lar mass versus NUV —r color and the BPT diagram 
(Baldwin et al. 1981). We use estimates of stellar mass 
and NUV—r from the New York University Value-Added 
Galaxy Catalog (NYU-VAGC 24 ; Blanton et al. 2005a) 

24 http://sdss.physics.nyu.edu/vagc/ 
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Fig. 1.— P-MaNGA galaxies are plotted as large red or blue symbols on the stellar mass versus NUV —r color plane (left panel) and 
the BPT (Baldwin et al. 1981) diagram (right panel). The different fields of the P-MaNGA observations are distinguished by the different 
symbols as indicated, while the number within the symbols indicates the IFU Bundle for each galaxy. Galaxies are divided into two 
classes: centrally quiescent (red solid symbols) and star-forming (blue open symbols) galaxies, defined as those with D n (4000)> 1.6 and 
D n (4000)< 1.6 in the central spaxel of the P-MaNGA datacube. Distributions of a volume-limited sample selected from the SDSS are 
shown as grey scale for comparison. The dashed horizontal lines in the left-hand panel are for NUV — r =4 and 5, which are commonly used 
in the literature to divide galaxies into blue-cloud, green-valley and red-sequence. The dashed and dotted lines in the right-hand panel are 
the AGN classification line of Kauffmann et al. (2003b) and Kewley et al. (2006). 
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Fig. 2.— P-MaNGA galaxies are plotted as large red or blue symbols on planes of 4000A-break vs. H<5^4 (left panel) and H(5^4 vs. 
log 10 EW(Ho:) (right panel). Shown in solid and dashed lines are solar metallicity models of continuous star formation histories and star 
formation bursts calculated using the stellar population synthesis code of Bruzual h Chariot (2003). Distributions of a volume-limited 
sample selected from the SDSS are shown as grey scale for comparison. Symbols/colors are the same as in the previous figure. The 
dashed vertical line in the left-hand panel indicates a constant D n (4000)=1.6, while the dashed horizontal line in the right-hand panel is 
for EW(Ha)=2A; these are the typical values adopted in the literature to divide a galaxy into quiescent and star-forming populations. 


and NASA Sloan Atlas (NSA) 25 . NSA is a catalog of 
images and parameters of local galaxies based on data 
from SDSS (York et al. 2000), GALEX (Martin et al. 
2005) and 2MASS (Skrutskie et al. 2006); for details 
see Blanton et al. (2005b,a, 2011). The stellar mass for 
each galaxy is estimated based on its spectroscopically- 
measured redshift and the five-band Petrosian magni¬ 
tudes from SDSS photometric data, as described in de¬ 
tail in Blanton & Roweis (2007). The NUV —r color is 

25 http://www.nsatlas.org 


defined by the integrated light in the NUV-band from 
GALEX and the Petrosian magnitude in the r-band 
from SDSS, both corrected for Galactic extinction. Mea¬ 
surements of the emission line ratios for the BPT dia¬ 
gram are taken from the MPA/JHU SDSS database 26 
(Brinchmann et al. 2004). 

Figure 2 displays the galaxies in the D„(4000) versus 
EW(H(5^) and EW(H<Ui) versus log 10 EW(Ho:) planes, 
where D„(4000) is the 4000A break in the optical spec- 

26 http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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trum as defined by Balogh et al. (1999), EW(H<5^) the 
Lick/IDS index of the Hd absorption line defined by 
Worthey & Ottaviani (1997), and log 10 EW(Ha) the log¬ 
arithm of Ha emission line equivalent width. We use 
measurements of these parameters from the MPA/JHU 
database, which were obtained from the SDSS 3"0- 
diameter spectra, thus probing the central 1-2 kpc for 
the P-MaNGA galaxies. 

In both these figures, we use squares and circles for 
galaxies in Fields 4 and 9 respectively, and label each 
galaxy by the BundlelD used in those fields. We divide 
these galaxies into two subsets according to the 4000A 
break in the central spaxel of their datacube (see § 2.3): 
centrally quiescent (CQ) with D n (4000)> 1.6 and cen¬ 
trally star-forming (CSF) with D„(4000)< 1.6. The two 
subsets are highlighted in the figures with red solid and 
blue open symbols respectively. We will use the same 
symbols and colors for subsequent figures unless other¬ 
wise stated. For comparison, we have selected a volume- 
limited sample of 21,328 galaxies with stellar mass above 
1O 9 M 0 and redshift in the range 0.01 < 2 < 0.03 from 
the NSA. Distributions of this sample are plotted in gray¬ 
scale maps in Figures 1 and 2. 

As can be seen from Figure 1, the P-MaNGA sam¬ 
ple (though small in size) spans a wide range in both 
the color-mass and BPT diagrams, similar to that of the 
general population. In the literature, galaxies are usu¬ 
ally classified into three types according to their NUV —r 
color: red-sequence (NUV —r > 5), green-valley (4 < 
NUV —r < 5) and blue-cloud (NUV—r < 4). Accord¬ 
ingly, half of the sample galaxies fall in the blue-cloud 
region, indicative of ongoing star formation and rela¬ 
tively rich cold gas content. For all the galaxies ex¬ 
cept ma002_9 and ma003_9, the color-based classification 
and the D„(4000)-based classification are consistent with 
each other, in the sense that the blue-cloud galaxies are 
centrally star forming, while the objects in the green 
valley and red sequence combine to form the centrally- 
quiescent galaxy class. This result indicates that green- 
valley and red-sequence galaxies share similar properties 
in their central regions, and that their different global 
color can be mainly attributed to the different stellar 
population and gas content in their outer parts. 

One of the two exceptions, ma003_9, is a strong Seyfert- 
type AGN according to its location in the BPT diagram, 
and is an outlier in the color-mass diagram due to its 
relatively blue color at its stellar mass. As shown be¬ 
low, this galaxy presents a red core and a blue outer disk 
in the optical image, such that its classification would 
be more reasonable if the measurement of NUV —r color 
were confined to its central region. The strong contra¬ 
diction between the central and global classification in 
this case highlights the importance of spatially-resolved 
measurements. The other exception, ma002_9, is a dusty, 
inclined spiral. It is a centrally star-forming galaxy ac¬ 
cording to its small value of D n (4000), but is classified as 
a green-valley galaxy based on the relatively red NUV —r 
color, which in this case is likely caused by dust. 

The D„(4000) and LWa indices are known to be indica¬ 
tors of the recent star formation history of galaxies, with 
D n (4000) being a good proxy for the luminosity-weighted 
stellar age, and Hd^ a sensitive tracer of the star forma¬ 
tion that occurred 0.1 1 Gyr ago (Bruzual & Chariot 


2003, hereafter BC03). As shown in Kauffmann et al. 
(2003b), the location of a galaxy in the plane of D„(4000) 
and H<5_4 provides a powerful diagnostic of whether the 
galaxy has experienced a continuous star formation his¬ 
tory or bursts over the past 1-2 Gyr. This property 
is illustrated in the left-hand panel of Figure 2, where 
we present the D„(4000) and indices for the BC03 
models of solar metallicity that follow exponentially de¬ 
clining star formation histories (SFR oc exp(—t/r)), ei¬ 
ther continuous star formation decline with long e-folding 
times (r > 5 x 10 8 yr, solid lines), or bursts of star for¬ 
mation with fairly short e-folding times (r < 5 x 10 8 yr, 
dashed lines). Different colors indicate the adopted val¬ 
ues of the characteristic timescale r. The right-hand 
panel presents the same sets of models on the HJ^ in¬ 
dex versus log 10 EW(Ha) plane. The Ha luminosity is 
computed from the output of the BC03 models. We 
convert Lyman continuum photons to Ha photons fol¬ 
lowing Hunter & Elmegreen (2004, see equations B2-B4 
in their appendix). We take the recombination coeffi¬ 
cients and Ha/H (3 ratios from Hummer & Storey (1987). 
For the LMC, solar, and super-solar metallicity mod¬ 
els we assume nebular temperatures of T e =15,000 K, 
10,000 K, 5000 K respectively. The P-MaNGA galaxies 
occupy roughly the same region as the general popula¬ 
tion in both panels, with the CSF galaxies located in the 
low-D„(4000) and high-Hd^ region and the CQ galaxies 
in the opposite corner. It is, however, notable that the 
P-MaNGA sample lacks galaxies of intermediate values 
in these indices, and consequently the current sample is 
somewhat biased to a more bimodal distribution on both 
planes compared to the full population of SDSS galaxies. 

In the rest of this paper we will use the D n (4000)- 
EW(H(5 j 4) and EW(H(5 J 4)-log 10 EW(Ha) planes to inves¬ 
tigate the maps and radial profiles of recent star for¬ 
mation histories of the P-MaNGA galaxies. In this 
case we use relative measurements of the three in¬ 
dices as indicators of recent star formation history, 
but we do not present quantitative measurements of 
mean stellar age, which depends on the sensitivity 
of these indices to stellar metallicity, element abun¬ 
dance and dust attenuation (e.g. Worthey & Ottaviani 
1997; Bruzual & Chariot 2003; Thomas et al. 2004; 
Thomas & Davies 2006; Sanchez et al. 2012a). In a par¬ 
allel paper (Wilkinson et al. 2015) we attempt to simul¬ 
taneously obtain the stellar age, metallicity and dust con¬ 
tent for the P-MaNGA galaxies, by fitting stellar popu¬ 
lation models to the full spectrum. In addition, in the 
present paper we assume that Ha emission is just a signa¬ 
ture of recent star formation. However, EW(Ha) lower 
than 6A was found to be more compatible with light 
from post-AGB stars (Papaderos et al. 2013; Sarzi et al. 
2013). One should keep these caveats in mind when inter¬ 
preting the observational results presented in this paper. 

2.3. Spectral fitting 

For each spectrum in the datacubes generated as de¬ 
scribed in §2.1, we have performed a decomposition of 
the emission-line component and the continuum plus 
absorption-line component (hereafter called the ‘stellar 
component’), using both the public spectral fitting code 
STARLIGHT (Cid Fernandes et al. 2004) and our own code 
described in Li et al. (2005). We have carefully masked 
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Fig. 3. — Example spectra of the P-MaNGA galaxies. The top two rows of panels are for a high-S/N spectrum with S/N ~ 20 at 5500A 
in the continuum, while the bottom two rows of panels are for a spectrum with S/N ~ 8 defined in the same way. In each case the upper 
panel shows the whole wavelength range, while the lower three panels show the three wavelength intervals with prominent emission lines. 
In each panel the black and red lines are the observed and best-fit spectra, while the blue line in the lower smaller panel is the difference 
between the two. The vertical, gray-shaded regions indicate the bandpasses of the emission-line region which is masked during the spectral 
fitting. 
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out the emission-line regions when doing the spectral 
fitting. Figure 3 displays the results of this procedure 
when applied to two representative spectra from the 
MaNGA datacubes, one with a high signal-to-noise ratio 
(S/N ~ 20) at 5500A in the continuum and one with a 
more modest S/N ~ 8. 

We measure the D„(4000) and indices from each 
best-fit spectrum, defining the indices in the same way 
as in the previous subsection. We have also measured 
the flux and equivalent width of the emission lines (both 
Balmer and forbidden) by fitting a Gaussian profile to 
these lines in the stellar component-subtracted spectrum. 
The emission-line parameters from STARLIGHT were com¬ 
pared to those obtained from our code. We have also 
performed the same procedure of spectral fitting and 
parameter measuring on Voronoi-binned versions of the 
datacubes to assess the effects of spatial averaging before 
fitting. All these analyses were found to return indistin¬ 
guishable results, and so in what follows we will only 
present the values we obtain when applying our code to 
the unbinned datacubes. 

We have also applied our code to the archival SDSS 
spectra of the centers of the P-MaNGA galaxies, and 
in Figure 4 we compare the resulting measurements of 
D„(4000), EW(Hdyi) and log 10 EW(Ha) to those ob¬ 
tained from the central spaxel of the P-MaNGA dat¬ 
acubes (plotted as blue squares in the figure). In the 
same figure, we also compare our measurements, from 
both SDSS (black crosses) and the central spaxel of P- 
MaNGA datacubes (red triangles) with the measure¬ 
ments taken from the MPA/JHU database. Overall, the 
values from different observations and pipelines agree 
well with each other. We have visually examined the 
best-fit spectrum for the few outliers, and find that our 
code provides a reasonably good fit to both SDSS and 
P-MaNGA spectra even in these cases. The relatively 
large discrepancy in one of the outliers (the blue square 
and black cross located in the left-bottom corner in the 
right-hand panel) is caused by the bad data pixels over 
the Ha band in the SDSS spectrum. 

We note that, when fitting the spectra and measuring 
the diagnostic parameters, we do not take into account 
the covariance between wavelengths, nor the covariance 
between pixels. Methods for tracking such covariances 
have been studied in depth (e.g. Sharp et al. 2015, for 
the SAMI survey). As pointed out in Wilkinson et al. 
(2015), when compared to P-MaNGA, MaNGA will pro¬ 
vide more accurately estimated spectral errors (thus ac¬ 
curate error propagation through to the final datacubes), 
allowing construction of covariances between pixels. We 
will apply a detailed treatment of covariance in MaNGA 
data in future work. 

3. RESULTS 

3.1. Global versus central parameters 

Before investigating the spatially-resolved properties 
of these galaxies’ star-formation histories, we take a first 
look at any spatial variations by comparing the central 
values of the diagnostic parameters to their integrated 
values across the field of view. To this end, we generated 
a global spectrum for each galaxy by stacking the entire 
datacube, with each spaxel at given wavelength being 
weighted by S/N 2 following Cappellari & Copin (2003), 


where S/N is the signal-to-noise ratio at 5500A in the 
continuum. Spaxels with S/N <3 were excluded from 
the stacking. We then measure the diagnostic parame¬ 
ters D„(4000), EW(Ma) and log 10 EW(Ha) in the same 
manner as above. Figure 5 compares the global measure¬ 
ments to the same parameters obtained from the central 
spaxel of the galaxies. For CSF galaxies, the global pa¬ 
rameters differ little from the central ones. In contrast, 
however, most of the CQ galaxies show significant differ¬ 
ence between central and global parameters, with smaller 
D„(4000) and larger EW(H<5 a) and log 10 EW(Ha) in the 
latter. This result implies that the outer regions of CQ 
galaxies contain younger stellar populations and have 
therefore experienced more recent star formation than 
their inner parts. The difference in EW(Ha) is relatively 
weak compared to the other two parameters, suggesting 
that the central-to-global variation may be mainly driven 
by the radial variation in recent star formation history, 
but not the current star formation rate. It may also be 
the case that Ha is not a reliable tracer of star forma¬ 
tion in galaxies dominated by older stellar populations 
(see § 4). 

A comparison of the global parameters with the cen¬ 
tral parameters on the diagnostic diagrams of recent star 
formation history is shown in Figure 6, where we plot 
the galaxies on the D„(4000)-EW(H<5 a) and EW(Hi5 j 4)- 
log 10 EW(Ha) planes, using the central and global pa¬ 
rameters for the upper and lower panels, respectively. 
Comparing the upper and lower panels, we identify two 
noticeable effects. First, the gap between the CSF and 
CQ galaxies as seen in the upper panels is largely filled 
in the lower panels, and this in-fill is made up exclusively 
of CQ galaxies due to the broad scatter in their global 
values of D„(4000) and EW^^)- Second, when glob¬ 
ally measured, both the CSF galaxies and the CQ galax¬ 
ies appear to deviate from the continuous star-formation 
models, moving by varying degrees toward the regime of 
recent bursts. 

3.2. 2D Maps and radial profiles 

In Figures 7 and 8 we present the 2D maps of the 
three diagnostic parameters, D„(4000), EW(Hd^) and 
log 10 EW(Ha), for all the 12 galaxies in the sample, with 
the CSF galaxies in one figure and the CQ galaxies in the 
other. In each case, the galaxies have been ordered such 
that the central values of D„(4000), as measured from 
the central spaxel of the datacubes, increases from left to 
right. The SDSS gri image is also shown for each galaxy. 
The effective radius, length and direction of major and 
minor axes are indicated in the images and maps. When 
generating the maps, we only employ spaxels where the 
S/N in the continuum at 5500A is greater than three. 

Consistent with what was found from the global mea¬ 
surements in the previous subsection, the overall impres¬ 
sion from the 2D maps is that all the galaxies show vari¬ 
ation within the probed radius in all three parameters, 
but to varying degrees and differing radial extents, with 
the variations more pronounced in the CQ galaxies. As 
noted above in §2.2, galaxy ma003_9 shows a blue outer 
disk in the optical image, which causes it to be classi¬ 
fied as a blue-cloud galaxy in the color-mass diagram, 
although its red core suggests that this is intrinsically a 
CQ galaxy. In agreement with this interpretation, the 
galaxy shows strong contrast between the central and 
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Fig. 4. — Measurements of D n (4000) (left panel), EW(H<5>i) (center panel) and log 10 EW(Ho') (right panel) for the P-MaNGA galaxies. 
Blue squares compare the measurements obtained by applying our code to the central spaxel of P-MaNGA datacubes and the SDSS spectra, 
which are compared in red triangles and black crosses to the measurements taken from the MPA/JHU database. For clarity errors are 
shown only on the blue squares, which compare the measurements obtained from the SDSS and P-MaNGA spectra. For given parameter 
its error is estimated by the 1 — a scatter between the three measurements. 
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Fig. 5.— Upper panels compare the global D n (4000), EW(H<5 A ) and log 1 QEW(Ho) of the P-MaNGA galaxies with the same parameters 
measured from the central spaxel of their datacubes. Lower panels show the difference between the global and central parameters as a 
function of the central one. The global parameters are measured from the average spectrum for each galaxy obtained by stacking the entire 
datacube, with each spaxel being weighted by S/N 2 . Symbols and colors are the same as in Figure 1. For given parameter, the error is 
given by the l-c scatter between the three different measurements as shown in Figure 4. 


outer parts in the maps of all the three diagnostic pa¬ 
rameters. Similar behavior is seen in galaxy ma008_4, 
the most massive galaxy in the P-MaNGA sample with 
M» = 2 x 10 11 Mq, which is also a face-on spiral with a 
relatively blue outer disk (although less prominent than 
in ma003_9) and strong radial variation in the diagnos¬ 
tic parameters. The galaxy is classified as a green-valley 
galaxy on the color-mass diagram (see Figure 1) due to 
its intermediate color, NUV— r= 4.4, which is unusually 
blue for its large stellar mass. 

The radial variations of the galaxies are shown more 
clearly in Figure 9, where we plot the radial profiles of 
the three parameters for all 12 galaxies. To save space 
we display two galaxies in every panel, with one CSF 
galaxy (blue symbols) and one CQ galaxy (red symbols). 
This paired format also highlights the highly different 
behavior of the two broad classes. In constructing these 
profiles, we have corrected for the effects of inclination 
on the radius for each spaxel in the datacubes, using the 
minor-to-major axis ratio from NSA, as determined from 
the SDSS r-band photometry. When estimating the ra¬ 


dial profile, we bin all the spaxels with the continuum 
S/N > 10 at 5500A into a set of radial intervals with a 
constant width of Alog 10 (R/R e ) = 0.2, according to the 
deprojected radius of the spaxel (r) and the effective ra¬ 
dius of the galaxy (Re). The value of the profile at given 
radius is then estimated by the median of the spaxels 
falling in the radial interval, and the error is given by 
the la scatter between the spaxels. 

Figure 9 clearly shows that CQ galaxies generally 
present significant radial gradients in all the parameters, 
in the sense that the outer part of the galaxies show 
weaker D n (4000) (thus a larger fraction of stellar pop¬ 
ulations with ages of < 1 Gyr), stronger Hd absorption 
(thus more recent star formation in the past a few Myr), 
and stronger Ha emission (thus likely stronger ongoing 
star formation). In contrast, CSF galaxies possess quite 
flat profiles in all cases, revealing little by way of gra¬ 
dient in star-formation properties. At least out to R e , 
the radial profiles are close to linear for all the galax¬ 
ies. For each galaxy we apply a linear fit to the ra- 
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Fig. 6. — Distribution of the P-MaNGA galaxies (symbols) on the plane of 4000A-break vs. index (left-hand panels), and the plane 
of H<5 a index vs. log 1 QEW(Ho), the equivalent width of Ha emission line in logarithm (right-hand panels). Results in the upper panels are 
for the central spaxel of the P-MaNGA datacubes. Results in the lower panels are based on the average spectrum, obtained from stacking 
the whole datacube for given galaxy. Shown in solid and dashed lines are solar metallicity models of continuous star formation histories and 
star formation bursts calculated using the stellar population synthesis code of Bruzual & Chariot (2003). Symbols for the MaNGA galaxies 
are the same as in the previous figure. For given parameter, the error is given by the 1-er scatter between the three different measurements 
as shown in Figure 4. 


dial profile within R e and show the best fit as a solid 
line in the figure. The slope of the best-fit line is in¬ 
dicated in each panel. Figure 10 compares the slope in¬ 
dex for the three parameters, a(D ra (4000)), a(EW(H<5yi)) 
and a(log 10 EW(Ha)). These plots confirm that the ra¬ 
dial gradients in the CSF galaxies are small, but also 
shows that the different parameters are well correlated 
across galaxies of both types. This correlation is partic¬ 
ularly strong for a(D„(4000)) and a(EW(H5 J 4)), which 
are beautifully (albeit non-linearly) correlated. However, 
given the small sample size, these apparent correlations 
should probably not be over-interpreted. 

Figure 11 compares the slopes of the radial profiles in 
a different way, where we plot all the profiles in the same 
panel for a given parameter, normalizing each profile by 
the inner-most radial bin and scaling the radii by the ef¬ 
fective radius. The differing behavior of the two classes of 
galaxies is striking. The CSF galaxies show only weak or 
no radial variation, while the CQ galaxies display steeper 
profiles, spanning a wide range in slope. It is, however, 
interesting to note that although the average properties 
of the two types differ dramatically, there is a degree of 
overlap between the classes in all panels, suggesting that 


a continuum of the same underlying physical processes 
has driven star formation in all these galaxies. 

As indicated from the 2D maps, the galaxies with 
a blue outer disk (ma003_9 and ma008_4) indeed show 
the strongest radial variation in all parameters. These 
two galaxies manifest their peculiarity in the profile of 
log 10 EW(Ha), which is flat (for ma008_4) or even neg¬ 
ative (for ma003_9) at R < 0AR e , before rapidly in¬ 
creasing at larger radii. When compared to these two 
galaxies, the other CQ galaxies present much shallower 
profiles (though still steeper than that of CSF galaxies), 
with log 10 EW(Ha) increasing smoothly at R < 0.6 R e 
and becoming flat at larger radii. The large values 
of log 10 EW(Ha) at large radii indicates that these two 
galaxies are strongly forming stars in their outer disks. 
This measurement is consistent with blue NUV—r color 
in these outer disks, which is known to be a sensitive 
indicator of the cold gas mass fraction in galaxies (e.g. 
Catinclla et al. 2010; Li et al. 2012a). 

The negative log 10 EW(Ha) profile at R < 0.4i? e seen 
in ma003_9 is also found in ma003_4, and may be caused 
by the presence of a bulge or AGN in the galactic cen¬ 
ter. Indeed, according to the BPT diagram (Figure 1), 












Gradients in recent star formation histories 


11 


ma004_4, R e =2.71 



D4000 






ma001_9, R e =2.90 






ma002_4, R e =3.90 




Fig. 7.— Panels from left to right show the SDSS image and maps of D n (4000), EW(H^) and log 10 EW(Ho:) for the centrally star¬ 
forming galaxies in P-MaNGA. From top to bottom panels the galaxies are sorted by increasing the D n (4000) measured from the central 
spaxel of the P-MaNGA datacube. The longer black line in each panel indicates the direction of the major axis, with the line length equal 
to twice the effective radius of the galaxy ( R e ). The shorter line perpendicular to the major axis indicates the minor axis, and the line 
length is set to b/a times R e , where b/a is the minor-to-major axis ratio determined from the r-band SDSS image. The P-MaNGA galaxy 
name, as well as the value of R e is indicated above the left-most panel, where the size of the SDSS images is set to 6R e x6R e . The hexagon 
in each panel indicates the field of view of the IFU bundle. The color scales in the D n (4000), EW(H5^) and log 10 EW(Ho;) maps are 
tailored to each galaxy for the sake of contrast. 
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ma003_9, R s =20.9 
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Fig. 8.— SDSS image (left-most panels) and maps of D n (4000), EW(H<5yt) and log;|nEW(Ha) for the centrally quiescent galaxies in 
P-MaNGA. The galaxies from top to bottom are ordered by increasing their central D n (4000). Symbols and lines are the same as in the 
previous figure. 
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Fig. 9.— Radial profiles of three diagnostic parameters measured for the P-MaNGA galaxies (from left to right): D n (4000), EW(H5 a) 
and log 10 EW(Ho:). Each panel shows results for two galaxies, a centrally star-forming galaxy in blue and a centrally quiescent galaxy 
in red, with their name (in form of Unit_Field) and global NUV-r indicated. From top to bottom panels the galaxies are ordered by 
increasing their central D n (4000). A linear fit to each radial profile within the effective radius (Re) is plotted as the solid line, with the 
slope index a indicated. 


the two galaxies are the only AGN candidates in the 
sample. It is likely that an AGN contributes to the Ha 
emission in the central region, thus enhancing the central 
log 10 EW(Ha) and leading to a dip in the profile as seen 
at R ~ 0.3I? e . Interestingly, apart from this feature, the 
radial profiles of the two probably AGN hosts are similar 
to those of the other CQ galaxies, suggesting that cur¬ 
rent, low levels of AGN activity may not have a strong 
effect on recent star formation activity, at least for these 
two examples. 

It is interesting to note that the three galaxies dis¬ 


cussed above, ma003_9, ma008_4 and ma003_4, share a 
number of common properties: they are all spiral galaxies 
with a quiescent center, at the massive end of our sample 
with M*~ lO n M0, and showing the largest gradients in 
D n (4000) indicative of young stellar populations formed 
within the past a few Gyr. They all share some properties 
with the population of “passive red spirals”, identified 
recently as an interesting set of possible transition ob¬ 
jects with disk-like morphologies but red colors, at both 
low-z (e.g. Wolf et al. 2005, 2009; Bamford et al. 2009; 
Masters et al. 2010) and high-z (e.g. Bundy et al. 2010). 
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Fig. 10.— Comparisons of the slope index, determined from the radial profile at radii within the effective radius for the three diagnostic 
parameters: D n (4000), EW(H5 J 4) and log 10 EW(Ho:). 





Fig. 11.— Radial profiles of three diagnostic parameters measured for the P-MaNGA galaxies (from left to right): D n (4000), EW(H^a) 
and log 10 EW(Ho:). For a given parameter and a given galaxy, each profile is normalized by the inner-most radial bin and the radii are 
scaled by the effective radius. Results for centrally star-forming and centrally quiescent galaxies are plotted in blue and red symbols/lines. 
The big symbols indicate the field and IFU unit of the P-MaNGA observation. 


Statistical analyses on a larger sample from MaNGA (or 
other IFU surveys) is needed in order to better under¬ 
stand the physical link between our galaxies and the pas¬ 
sive/quiescent spiral population. 

The 2D maps in Figures 7 and 8, as well as the ra¬ 
dial profiles in Figures 9 and 11, are in broad agreement 
with maps/profiles of stellar population age estimated 
from both CALIFA (Gonzalez Delgado et al. 2014) and 
P-MaNGA (Wilkinson et al. 2015), in the sense that 
most galaxies in the nearby universe show flat to nega¬ 
tive gradients in stellar age, although in those studies the 
detailed behaviors of the age gradients depend on galaxy 
mass or morphology. In this work we only consider the 
dependence of our measurements on the D„(4000)of the 
galactic center, which is an indicator of the prominence of 
stellar populations younger than 1-2 Gyr, rather than the 
mean stellar age as investigated in those studies. Despite 
these differences, the IFU-based results obtained so far 
are all consistent with the “inside-out” growth of galac¬ 
tic disks (White & Frenk 1991; Mo et al. 1998), a picture 
that was supported previously by both numerical simu¬ 
lations (e.g. Brook et al. 2006) and abundant studies of 
broadband color gradients pioneered more than 50 years 
ago by Tifft (1963). 


3.3. Diagnostic diagrams of recent star formation 
history 

We now return to the comparison of diagnostic dia¬ 
grams of recent star formation history, based on the ra¬ 
dial profiles obtained in the previous subsection. Fig¬ 
ure 12 presents the D„(4000) vs. EW(H^a) (left panels) 
and EW(H<5 a) vs. log 10 EW(Ha) (right panels) planes, 
with each panel once again contrasting two galaxies, one 
CSF system (blue symbols) and one CQ system (red sym¬ 
bols). Panels are ordered such that D„(4000) increases 
from top to bottom. For each galaxy, the radial profiles 
within the effective radius R e as shown in the previous 
figure are plotted here, with larger data points corre¬ 
sponding to larger radii. The set of BC03 models in 
Figure 6 are also repeated here as the solid and dashed 
lines, but for clarity these lines are not colored; they pro¬ 
vide a reference grid to register the relative position of 
the different radial bins on the diagnostic diagrams. 

Both the weak radial variations for CSF galaxies and 
the strong variations for CQ galaxies are readily ap¬ 
parent in this figure. Comparing the location of the 
radial bins on the two diagrams, as well as their rela¬ 
tive location with respect to the model curves, presents 

















8 

6 

4 

2 

0 

§ 

8 

6 

4 

2 

0 

§ 

8 

6 

4 

2 

0 

§ 

8 

6 

4 

2 

0 

§ 

8 

6 

4 

2 

0 

§ 

8 

6 

4 

2 

0 

2 

1 

2 , 


Gradients in recent star formation histories 


15 



1.9 2.0 



EW(H<5 a ) 


— Radial variation of the P-MaNGA galaxies in the plane of D n (4000) vs. EW(HSa) (left panels), and the plane of EW(Hi5a) 
W(Hct) (right panels). Each panel shows results for two galaxies: a centrally star-forming galaxy (blue symbols) and a centrally 
jalaxy (red symbols), which are defined to have D n (4000)< 1.6 and On(4000)> 1.6 in their center. Larger data points correspond 
idius. Overplotted in solid and dashed lines are the same set of stellar synthesis models as in previous figures. Galaxy name and 
re indicated in each panel. 
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a comprehensive picture of how and where star forma¬ 
tion cessation is occurring within the galaxies. The di¬ 
agrams demonstrate that all the CQ galaxies behave in 
the same manner: D n (4000) decreases, while EW(H<5 a) 
and log 10 EW(Ha) increase, as one moves outward in the 
galaxies. This consistency is even clearer in Figure 13 
where we overplot the radial profiles of all 12 galaxies in 
a single set of diagnostic diagrams. The profiles are plot¬ 
ted with cyan triangles for CSF galaxies, and magenta 
triangles for CQ galaxies. We highlight the two peculiar 
galaxies discussed above, (ma003_9 and ma008_4), which 
show the strongest radial gradients and unusual shape 
in the profile of log 10 EW(Ho:), by identifying them with 
yellow triangles. We also pick out the green-valley galax¬ 
ies (4 <NUV—r< 5) as green triangles. For comparison, 
the central spaxels are plotted with blue open symbols 
and red solid symbols for the two classes of galaxies. The 
CQ galaxies indeed follow a tight track on the D„(4000) 
versus EW(HcM) plane. 

These spatially-resolved data, reaching out to R e for 
all 12 galaxies, also show how the clear distinction be¬ 
tween the two classes of galaxies becomes blurred when 
one also considers their outer parts. Figure 13 shows 
that the gaps in the parameter space between the re¬ 
gions occupied by the centers of the galaxies are popu¬ 
lated by these same galaxies’ outer parts. This contin¬ 
uous sequence forms a tight relation in the D„(4000)- 
EW(ffiyi) plane, closely following the continuous star 
formation models. In detail, though, the intermediate 
region on the diagrams is mainly occupied by the two 
galaxies (ma003_9 and ma008_4) that have strong radial 
variations and unusual profiles in log 10 EW(Ha). In or¬ 
der to see this dominance more clearly, we have repeated 
the same panels of Figure 13, but excluding these two 
systems, and we find the gap between the galaxy types 
once again starts to emerge. It is clear, however, that the 
CQ galaxies do still extend substantially into the inter¬ 
mediate region, and that the effect appears to be stronger 
for green-valley galaxies than for the other objects in the 
CQ galaxy class. 

3.4. Notes on individual galaxies 

ma008_4: This system, CGCG 184-033, is the central 
galaxy of a group in the catalog of Yang et al. 
(2007). The system comprises 8 member galax¬ 
ies with spectroscopically-measured redshift from 
SDSS, including the central galaxy itself. This is 
the most massive galaxy in the sample with a stel¬ 
lar mass of M* = 2 x 1O 11 M 0 . The SDSS spec¬ 
troscopy and the central spaxel of the P-MaNGA 
datacube reveal that the central region of this 
galaxy shows relatively strong absorption in H<5 for 
its EW(Ha), implying that a starburst event has 
happened within the galaxy in the past 1 - 2 Gyr. 
Considering the fact that the galaxy has a very 
close companion, a strongly star-forming spiral at 
the same redshift (z=0.042), the post-starburst fea¬ 
ture may be the product of a tidal interaction with 
the companion, which caused the central part of 
the galaxy to have ceased star formation recently. 
The P-MaNGA data also reveal strong radial varia¬ 
tion in all the diagnostic parameters in this galaxy: 
D„(4000), EW(H(5 a) and EW(Ho:). In particular 
the profile of log 10 EW(Ha) shows a very strange 


shape when compared to the other CQ galaxies. 

ma003_9: This galaxy is a spiral with a red core but 
very blue outer disk in the optical image (see Fig¬ 
ure 8). It is classified as a blue-cloud galaxy be¬ 
cause of its blue NUV —r color, but it is clas¬ 
sified as a centrally quiescent galaxy due to its 
large central D„(4000). The red core with old 
stellar populations and a blue global color seems 
to suggest that the galaxy might be a ‘Bluedisk’ 
galaxy, a class of system with an unusually high 
fraction of Hi gas mass, as recently studied in de¬ 
tail by Wang et al. (2013). However, the Hi gas 
mass is log lo (MHi/M 0 ) = 9.8 according to the Hi 
global profile from Huchtmeier & Richter (1989), 
homogenized in Paturel et al. (2003). This is com¬ 
parable to the value expected from the NUV —r 
and stellar surface mass density for a normal spi¬ 
ral (Catinella et al. 2010; Li et al. 2012b). The 
BPT diagram shown in Figure 1 reveals strong 
nuclear activity in the galactic center. A ~ 
12 mjy source coincident with this galaxy is de¬ 
tected by the NVSS survey (Condon et al. 1998), 
but no source is detected by the FIRST sur¬ 
vey (Becker et al. 1995), suggesting that the radio 
emission from this galaxy is variable. This is con¬ 
sistent with an AGN origin, albeit at a low luminos¬ 
ity (Ti.4GHz ~ 10 21 - 6 W/Hz). See the related pa¬ 
per by Belfiore et al. (2015) for a spatially-resolved 
BPT diagram of this galaxy, also based on the P- 
MaNGA data. 

ma003_4: This system is one of the two AGNs in the 
sample (the other one is ma003_9 above). The AGN 
component is likely responsible for the strange pro¬ 
files of the log 10 EW(Ha) diagnostic in these sys¬ 
tems, as discussed above. 

ma006_4: With a stellar mass of 6.3 x 1O 9 M 0 , this 
galaxy is the only low-mass CQ galaxy in the P- 
MaNGA sample. It has unusually low EW(Hq) 
for its EW(H5 J 4)(see Figure 12). This might be 
a metallicity effect, as the galaxy is much closer 
to the continuous models for a lower metallicity of 
Z=0.008. 

4. SUMMARY AND DISCUSSION 

Using datacubes for 12 galaxies produced during the 
MaNGA prototype (P-MaNGA) observations, we have 
obtained maps and radial profiles of the D„(4000), 
EW(H5^) and EW(Ha) spectral diagnostics. These pa¬ 
rameters were derived for each spaxel by decomposing 
the spectrum into an emission-line component and a con¬ 
tinuum plus absorption-line component, and measuring 
the relevant indices from these separate components. We 
study the spatially-resolved recent star formation histo¬ 
ries of these galaxies by displaying their radial profiles 
on the D„(4000) versus EW(H<5 a) and EW(H<5 a) ver¬ 
sus log 10 EW(Ha) planes. We have classified the galax¬ 
ies as either ‘centrally quiescent’ (CQ) or ‘centrally star¬ 
forming’ (CSF) according to the value of D„(4000) in the 
central spaxel, and we compare the 2D maps and radial 
profiles for the two classes. CQ galaxies present signifi¬ 
cant radial gradients in all the three parameters, in the 
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Fig. 13.— Same as Figure 6, except that the radial profiles of the P-MaNGA galaxies are additionally plotted, with magenta (light cyan) 
triangles for centrally quiescent (star-forming) galaxies. The profile for a given galaxy is connected with a solid line and the symbols are 
sized by radius, with larger triangles for larger radii. Plotted in yellow triangles are the two galaxies (ma003_9 and ma008_4), which show 
the strongest radial gradients and unusual shape in the profile of log 10 EW(Ho:), while green triangles highlight the green-valley galaxies 
with 4 <NUV—r< 5. The individual spaxels of the P-MaNGA galaxies are plotted in the background as grey scale. 


sense that D„(4000) decreases, while both EW(H^ J 4 ) and 
EW(Ha) increase as one goes from the galactic center to¬ 
ward the outer region. In contrast, CSF galaxies show 
very weak or no radial variations. The outer parts of the 
galaxies show greater scatter on the diagnostic diagrams 
than their central regions. In particular, the separation 
between CQ and CSF galaxies in these diagnostic planes 
is largely filled in by the outer parts of galaxies whose 
global colors place them in the green valley. 

The three spectral diagnostics are known to be sen¬ 
sitive indicators of stellar populations of different ages. 
Spatially resolving a galaxy on the D I j(4000)-EW(H(5 J 4)- 
EW(Ha) diagrams probes the radial variations in recent 
star formation histories, thus revealing the way in which 
the galaxy grows and its star formation ceases. When 
spatially resolved, the P-MaNGA galaxies are found to 
form a tight sequence on the D n (4000) versus EW(ffix) 
plane, although the distribution of their central regions 
is strongly bimodal (Figure 13). This sequence closely 
follows the continuous star formation locus predicted by 
current stellar population models, and covers a very sim¬ 
ilar area to the distribution of the large sample of galac¬ 
tic centers from SDSS (Figure 2). This relation suggests 
that, at least for the systems being studied here, galaxy 
growth has been a smooth process, and is likely regulated 
by a common set of physical drivers. This result is consis¬ 
tent with the ‘inside-out’ picture of galaxy growth, where 
the stellar mass assembly starts in the galactic center and 
gradually extends to the outer regions (White & Frenk 
1991; Mo et al. 1998; Brook et al. 2006). In this picture, 
the shutdown of star formation also first occurs in the 
central region and slowly propagates out to ever larger 
radii. However, this result should not be overemphasized 
given the small size of our sample. The tight sequence 
on the D„(4000)-EW(H<5 j 4) diagram might be just a re¬ 
sult of small number statistic. It could also be the case 
that there are a variety of physical processes behind this 
sequence, which can be discriminated only when a larger 


sample becomes available. 

Radially resolved stellar populations have been studied 
previously to some extent, based on multi-wavelength 
broadband photometry at both low-z (e.g. de Jong 1996; 
Bell & de Jong 2000; Taylor 2005; Munoz-Mateos et al. 
2007; Zibetti et al. 2009; Roche et al. 2010; Suh et al. 
2010; Tortora et al. 2010; Gonzalez-Perez et al. 2011; 
Tortora et al. 2011; Kauffmann 2015) and high- 
z (e.g. Abraham et al. 1999; Azzollini et al. 2009; 
Szomoru et al. 2012; Wuyts et al. 2012; Szomoru et al. 
2013; Wuyts et al. 2013), as well as long-slit spec¬ 
troscopy for nearby galaxies (e.g. Moran et al. 2012; 
Huang et al. 2013b). Most spirals are found to have 
stellar population gradients, with the inner regions being 
older and more metal-rich than their outer regions, sup¬ 
porting the inside-out scenario of disk galaxy formation. 
The majority of the early-type galaxies show flat color 
profiles, with only a small fraction being “blue-core” 
systems with negative color gradients (i.e. bluer outer 
parts). The star formation history of galactic disks is 
found to be primarily driven by its structural properties, 
as usually quantified by the surface stellar mass density 
(/x*, Bell & de Jong 2000) and stellar light concentration 
(Gonzalez-Perez et al. 2011), while the total stellar 
mass is a secondary parameter, correlated more closely 
with metallicity than with age. For massive galaxies, 
recent short-lived episodes or bursts of star formation 
occurring in their outer regions are required to interpret 
the radially resolved spectral indices as measured from 
long-slit spectroscopy (Huang et al. 2013b), a result 
that is consistent with the existence of many early-type 
galaxies in the nearby universe which display extended 
star-forming discs, as originally discovered by ultraviolet 
imaging (e.g. Yi et al. 2005; Kauffmann et al. 2007; 
Fang et al. 2012). The overall radial dependence of the 
diagnostic parameters presented in our work is appar¬ 
ently consistent with these previous findings, although 
the small sample size doesn’t allow the correlation of 
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the radial profiles with galaxy morphology, stellar mass 
and structural properties to be examined in a statistical 
manner. 

More recent studies using data from CALIFA have 
provided the most comprehensive results so far regard¬ 
ing the radial variations of the stellar population and 
star formation history of nearby galaxies (Perez et al. 
2013; Gonzalez Delgado et al. 2014; Sanchez et al. 2014; 
Sanchez-Blazquez et al. 2014). Perez et al. (2013) ana¬ 
lyzed the first 105 galaxies of the CALIFA survey, quan¬ 
tifying their spatially resolved history of stellar mass as¬ 
sembly and demonstrating how massive galaxies grow 
their stellar mass inside-out. Oxygen abundances esti¬ 
mated from Hu regions in 306 CALIFA galaxies present 
a common negative gradient in the oxygen abundance 
between 0.3 and 2 times R e , independent of morphol¬ 
ogy, the presence of bars, absolute magnitude or mass, 
provided that the analysis is limited to those galaxies 
without clear evidence of an interaction (Sanchez et al. 
2014). This result also supports the inside-out picture 
for the growth of galactic disks, and is consistent with 
the stellar age and metallicity gradients obtained by 
Sanchez-Blazquez et al. (2014) based on CALIFA data 
for a smaller sample of 62 nearly face-on spirals. 

Using a sample of 107 galaxies from CALIFA, 
Gonzalez Delgado et al. (2014) performed a detailed 
study of the radial structure of the stellar mass surface 
density (/i*) and stellar population age, as well as their 
dependence on the total stellar mass and morphology. 
In support of the inside-out formation picture, nega¬ 
tive gradients of the stellar age are present in most of 
their galaxies. The behaviors of the D„(4000) profiles in 
the P-MaNGA galaxies (see Fig. 11) are in good agree¬ 
ment with this result, although the D„(4000) is not a 
direct measure of stellar age, but rather an indicator 
of the prominence of young stellar populations. Con¬ 
firming the previous broadband photometry-based find¬ 
ings, the CALIFA data reveals that /i» is more important 
than M* in driving the star formation history in disks, 
while M* plays a more fundamental role in spheroidal 
systems. The CALIFA results, together with the results 
from broadband photometry and long-slit spectroscopy, 
strongly suggest that galaxy growth and death are driven 
by multiple parameters (at least the local surface den¬ 
sity, total stellar mass and morphology), and that further 
studies of much larger samples are needed in order to iso¬ 
late the role of each parameter, and the physical process 
behind it. The increasing MaNGA sample is expected to 
provide much better results in the near future. 

Another interesting result from our work is that, when 
spatially resolved, the distribution of galaxies on the di¬ 
agnostic diagrams is no longer bimodal, and the sepa¬ 
ration between the CQ and CSF populations is filled in 
by those galaxies that are classified as ‘green-valley’ ac¬ 
cording to their optical/UV color. In a recent comple¬ 
mentary study, Mendel et al. (2013) selected from SDSS 
a large sample of ‘recently quenched’ galaxies that have 
a significant fraction of young stars within the central 3" 
region, and found them to populate the region of green- 
valley galaxies in D n (4000) and optical colors, so their 
central parts seem to be at the same evolutionary stage 
as the outer parts of the P-MaNGA galaxies. The two 
studies combine to suggest that the green valley does, 
indeed, pick out a transition phase between star-forming 


and passive populations, but that different parts of dif¬ 
ferent galaxies enter the green valley at different times. 

There are a few caveats to this analysis that we 
will address in future work. We have used relative 
measurements of D ra (4000), EW(H<5.4) and EW(Ha) 
as indicators of recent star formation history. We 
have not presented quantitative measurements of mean 
stellar ages due to the sensitivity of these indices to 
stellar metallicity, element abundance and dust atten¬ 
uation (Worthey & Ottaviani 1997; Bruzual & Chariot 
2003; Thomas et al. 2004; Thomas & Davies 2006; 
Sanchez et al. 2012a). Throughout this paper we use so¬ 
lar metallicity models for comparison. We have com¬ 
pared the data to models with three different metal- 
licities (Z=0.4, 1, 2.5Z 0 ), finding that the dependence 
on metallicity is very small for younger stellar popu¬ 
lations, but become more important at older ages. In 
Wilkinson et al. (2015) we explore the possibility of si¬ 
multaneously deriving the age and metallicity of P- 
MaNGA galaxies using stellar population model fits 
to the full spectrum. In the present paper, we use 
the Ha emission line to characterize the instantaneous 
star formation rate. However, for galaxies with dom¬ 
inant old stellar populations, OB stars may not be the 
only source of ionizing photons (e.g. Cid Fernandes et al. 
2010; Sanchez et al. 2014). Indeed, EW(Ha) lower 
than 6A was found to be more compatible with light 
from post-AGB stars (Papaderos et al. 2013; Sarzi et al. 
2013), and it would be interesting to evaluate the con¬ 
tribution of such stars to the emission-line properties of 
the spectra once a larger sample is available. 

In a recent SDSS-based study Kauffmann (2014) has 
nicely demonstrated that a combination of D„(4000), 
EW(H^) and the specific star formation rate (SFR/M*) 
provides stringent constraints on the recent star forma¬ 
tion histories of the central region of galaxies in the stellar 
mass range 10 s — 1O 1O M 0 . In particular, this methodol¬ 
ogy is powerful in separating out the galaxies with ongo¬ 
ing or past starbursts from those with a continuous star 
formation history. By applying this new methodology 
to the spatially resolved spectroscopy for a large sam¬ 
ple from MaNGA, we may well expect to better quantify 
the fractions of the mass produced by ongoing and recent 
bursts, as well as their radial variation across the whole 
galaxy, thus leading to more reliable constraints on the 
maps and radial profiles of the recent star formation his¬ 
tory. 

Previous studies have proposed a variety of galaxy 
properties as indicators of the termination of star for¬ 
mation, including both color indices (e.g. NUV—r) and 
structural parameters like the central surface mass den¬ 
sity and stellar velocity dispersion (e.g. Fang et al. 2013), 
the effective radius R e and the ratio of stellar mass rel¬ 
ative to R e , R e 2 , or .Re 1 ' 5 (e.g. Omand et al. 2014), and 
Sersic profile index n (e.g. Bell 2008). Since different 
indicators may trace different physical processes, it is 
important to ask which galaxy property is most closely 
linked to the end of star formation. Again, once larger 
samples of galaxies become available, it will be very in¬ 
teresting to correlate the kind of spectral diagnostic gra¬ 
dient study presented here with these other parameters, 
to see which is the trigger, and even how the trigger 
varies between galaxies of different types or in different 
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environments. 

Despite these limitations, it is heartening to see how 
much can be learned about the cessation of star forma¬ 
tion from even a sample of 12 galaxies, thanks to the re¬ 
solved spectroscopic observations that P-MaNGA IFUs 
have provided. The full MaNGA survey began on 2014 
July 1, offering new opportunities to extend these kinds 
of studies, as observations of its 10,000 galaxies continue. 
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